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ABSTRACT 

A QCD relativistic potential model is employed to compute the decay rate 
and the photon spectrum of the process B~ — > pTv^. The result B(B~ — > 
fi^u^) ~ 1 x 10 -6 confirms the enhancement of this decay channel with 
respect to the purely leptonic mode, and supports the proposal of using this 
process to access relevant hadronic quantities such as the 5-meson leptonic 
decay constant and the CKM matrix element V u f,. 
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Noticeable theoretical attention has been recently given to the weak radiative decay 

B~ -> yrvfi . (1) 

The reason is in the peculiar role of this decay mode for the understanding of the dynamics 
of the annihilation processes occuring in heavy mesons [0], |2|, [3], |J. Moreover, it has been 
observed that ([!]) can be studied to obtain indications on the value of the 5-meson 
leptonic constant /# using a decay channel which differs from the purely leptonic modes 
B~ — > i~ i>g, and is not hampered by the limitations affecting those latter processes. Such 
difficulties mainly consist in low decay rates (using V u b = 3 x 1CT 3 , = 200 MeV and 
t b - = 1.646 ± 0.063 ps f§ one predicts B{B~ -> e~z> e ) ~ 6.6 10~ 12 and B(B~ -> // _ z/ M ) ~ 
2.8 10 ) or in reconstruction problems for .B — > t~v t . 

In ref. [Q heavy quark symmetry and experimental data on D*° — > D°7 have been 
exploited to study the dependence of B(B~ — > fjrv fi 'y) on the heavy meson decay constant 
F/y/mb, which is the common value of fs and fs* (modulo logarithmic factors) in the 
limit rrif, — > oo. The analysis is based on the dominance of polar diagrams contributing 
to the process B~ — > fi~v fl / -f, the pole being either the vector meson B* or the positive 
parity J p = 1 + state B[ (see ref. [Q for further details). According to the analysis in ||, 
in correspondence to the expected range of values of F: F ~ 0.35 GeV 3 ^ 2 , the branching 
ratio B(B~ — > jjTv^) should be (9(10 -6 ), which represents an enhancement with respect 
to the purely leptonic mode. 

In order to give further arguments in support of that analysis, we want to consider 
the process ([!]) in a different context. More precisely, whereas in we have studied 
the feasibility of extracting f% — F / ^frrTb from future experimental data, in this letter we 
study the decay (|]) within a well defined theoretical model in order to have an independent 
estimate of the decay rate. 

We employ a relativistic constituent quark model already used to study several aspects 
of the 5-meson phenomenology [8|, |J . Within this model the mesons are represented as 
bound states of valence quarks and antiquarks interacting via a QCD inspired istantaneous 
potential with a linear dependence at large distances, to account for confinement, and 
a modified coulombic behaviour at short distances to include the asymptotic freedom 
property of QCD. We adopt the interpolating form between such asymptotic dependences 
3 Present bounds are: B(B~ -> e~D e ) < 1.5 1(T 5 , B(B~ -> ^v^) < 2.1 1(T 5 § . 
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provided by the Richardson potential [TIJ Q. In the rest frame, the state describing a B a 
meson is represented as: 

l^>=*E%E% / dk 1 MhP(ki,r,a)d{(-k 1 ,s,P)\0>, (2) 

a /3 V O rs V ^ J 

where a and /? are colour indices, r and s are spin indices, &t and d\ are creation operators 
of the quark b and the antiquark q a , carrying momenta k\ and —ki respectively. The B- 
meson wave function ipB(ki) satisfies a wave equation with relativistic kinematics (Salpeter 
equation) JTTJ] taking the form (in the meson rest frame): 

{ Jkf + m\ + yjk\ + ml - M Ba }Mh) + Jdk ! 1 V(k l ,k' 1 )Mk[) = Q (3) 

where V(ki, k[) is the interaction potential in the momentum space and ips is covariantly 
normalized: 

J dhlMh)] 2 = 2 M Ba . (4) 

Solving eq.(^|) by numerical methods one obtains the wave function; we shall present 
this result later on. It can be mentioned that by this model a number of predictions 
have been derived; for example, the heavy meson spectrum, leptonic constants 0, ||, 
semileptonic form factors and strong decay constants 0. 

In the framework of the relativistic QCD potential model the wave function ipg is the 
main dynamical quantity governing the decay ([I]). As a matter of fact, the amplitude of 
the process B~(p) — > fi~(pi) v^ip?) j(k, e) can be written as 

A(B~ -> p-vpi) = ^V ub (L 1 * ■ n M ) , (5) 

where Gp is the Fermi constant, V u b is the CKM matrix element involved in the decay, 
L M = /i(p 1 )7 M (l — 7 5 )z/(p 2 ) is the weak leptonic current, and 11^ is defined by 11^ = U^ v e* u 
(e is the photon polarization vector); 11^ represents the correlator 

rV = i J d'xe^ < 0|T[J M (x)K(0)p(p) > . (6) 

In eq.(||) q is q = p\ + p 2 , J^{x) = u{x)^^{l — 7 5 )6(x) is the weak hadronic current and 

V u (0) = |e u(0)j u u(0) — |e 5(0)7„6(0) is the electromagnetic (e.m.) current. The two 

4 A smearing of the Richardson potential at short distances has also been introduced to take into 
account the effects of the relativistic kinematics; see ref.fl for the explicit form of the potential. 



pieces in the e.m. current correspond to the coupling to the light quark and to the heavy 
quark, respectively. The corresponding contributions to IL^ will be referred to as IFL, 

, , , . . ^ , , , , , , . . nnn tains t hp licrht rmar 

[IVl ^j/ivuvu "» "to"-" J-")"' "^j, 



and 11^ , depicted in figs. la,b. Ilf, contains the light quark propagator: 



e i£-x 

S u(x, 0) = / (27r)4g 2_ m 2 ^ + ™«) ( 7 ) 



while n^j, contains the analogous b quark propagator. 



The calculation of the time-ordered product appearing in (Q) gives: 

The analogous expression for 11^ can be obtained by the replacements: | — ~, m M «-> m fe , 
J^, <->• J* . The operators and J*,, which depend on the integration variable x, are 
written as J* = u(0)T^ u b(x) and = u(x)T^ u b(0), with the T matrices given by 
— lu(/+ ^«)7m( 1 _ 75) and = 7^(1 - 7 5 )(/+ m&^j,. In the constituent quark 
model such operators can be expressed in terms of quark operators; for example one has: 

J lv = £ £ <W 

a/3 rs 



r d 3 q 1 d 3 q 2 


r m b m u 


1/2 


1 (2vr) 3 


^„(gi)^(g 2 )- 





[IV 



[u u (q\, r)bl(q\,r, a) + v u (q\, r)d u (q\, r, a)] Y\ 
[u b (q 2 , s)b b (q 2 , s, (3)e-^ x + v b (q 2 , s)d\{q 2l s, (3)e^ x ] : (9) 



where E q (q) = ^jq 2 + vn? q and u q (v q ) are quark (antiquark) spinors. 

By exploiting anticommutation relations among annihilation and creation operators, 
we obtain, in the B meson rest frame: 



=*-7f£ / i^sMh) ^ (io) 



d 3 h 



m b m u 



1/2 



VE^J (2vr)3 r ^ E b (ki)E u (ki) 
r J 7^ 2 -^) +^]7^(l-75) 7^(l-75)[(^-^)+m & ]7 y | -. 
I (ft - <?) - (9i - Q) ~ m b J 



where q± = (E u ,—ki) and q 2 = (E b ,ki). We recognize in the two factors in the curly 
brakets the contributions of and 11^,, respectively. 

Since U^ u only depends on two vectors, the B meson momentum p^ and the photon 
momentum it can be written in terms of six independent Lorentz structures: 

IV = a Pf,p u + (3 k^ky + ( k^p v + 5 p^k u + £ g^ + irj e^ upa p p k a . (11) 
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By gauge invariance one has a = and £ = — p ■ k (; moreover, after saturation by e* u 
one gets: 

n M = n^e*" = [C (k^p u - p ■ kg^) + 1 7] e^pPk*) , (12) 

i.e. only the terms proportional to rj and ( survive: they are the vector and the axial 
vector contribution, respectively. It is convenient to compute eq.(^) in the B rest frame 
p = (M B , 0), with k = (k°, 0, 0, A; ); the result reads: 

nn nis 



C 



M B k° V i M B k° ' 

At this point, it is straightforward to calculate the rate of the decay process ([!]); one 
obtains (m M ~ 0): 



(13) 



T(B- -> irw) = £ lB/2 dk°k°(M B - 2fc°)[|n 11 | 2 + |n 12 | 



(14) 



where: 

nn 



i e 



dcosO / |A;i|d|A;i| u B (\ki\) 
Jo 



E b E u (E b + m b )(E u + m u )- 
(M B - k°)[(E b + m b )(E u + m u ) - \h\ 2 } + \k l \ 2 cos 2 6(E u +m u - E b - m b ) 



AV3tt J-i 
2 r 

7 

+ \ki\k°cos9(E u + m u — E b — m b ) 
- (E b + m b )(E u + m u )(E b + m u ) + (E b - m u )\ki\ 2 
1 r 
9 

\ki\k°cos9(E u + m u — E b — m b ) 



1/2 



(15) 



+ 



(M B - A; ) pi | 2 - (E b + m b )(E u + m u )} + \k l \ 2 cos 2 6(E u + m u - E b - m b ) 



+ (E b + m b ) (E u + m u ) (E u + m b ) - (E u - m b ) \ k\ 



and 



n 



12 



{ 



'" J /* |^1 1 max —t—t —> 

- , dcosO I \ki\d\ki\ u B (\ki\) 

4V37T J -i Jo 



1 



E b E u (E b + m b )(E u + m u ) 



1/2 



k°[(E b + m b )(E u + m u ) - \k x \ 2 } + \ki\cos9(m b - m u )(E u + m u + E b + m b ) 

2 1 

+ |^i|cose(M B -A: )(K + m u -^ 6 -m 6 )}(---) . (16) 



In eqs. ( |15| , |T6D the quantities / and g are defined as 

f = m 2 b + M 2 B - 2M B k° - 2M B E b + 2E b k° - 2A;°|£i|cos0 - m 2 (17) 



g = m 2 u + M 2 B - 2M B k° - 2M B E U + 2E u k° + 2k \h\cos9 - m\ , (18) 
whereas the S-wave reduced function u is related to ipB'- 

u B {\k l \) = ^-i, B {k 1 ) ■ (19) 

we plot in fig.2 the function u B obtained as a solution of the wave equation (|3]). 

Let us now point out that in computing the diagrams in figs.la,b and, therefore, in 
eqs. ([l5|) , (|16D , we have so far imposed 4- momentum conservation for the physical parti- 
cles B, fi, v and 7 in the process ([!]). On the other hand, energy conservation has to be 



imposed also at quark level since, otherwise, (|15| ) and (16) would present spurious kine- 



matical singularities. In order to deal with this problem we follow the approach originally 



proposed within the ACCMM model [12| for the decay b — > u £ i>£. One assumes that the 
spectator quark has a definite mass, while the active quark has a "running" mass, defined 
consistently with the energy conservation: 

E b + E U = M B . (20) 

Therefore, as in the case of the ACCMM model, the running mass of the active b quark 
can be defined by: 



mKh) = M 2 B + m 2 u - 2M B ^k\ + ml . (21) 

Moreover, by requiring that the right hand side of eq. (|21|) is positive, an upper bound 
on the quark momentum \k\\: can be obtained 

, r , Ml- ml , s 

N < — • 22 

1 1 - 2M B 

Notice that the masses of the light constituent quarks, as obtained by fits to the meson 
spectrum, are m u = m^ = 38 MeV. 

The contribution of the two physical processes when the photon couples to the light 
or to the heavy quark is still recognizable in eqs. (pT5|) - ([16D , since the quantities / and 
g come from the light and heavy quark propagators respectively. It turns out that the 
contribution of the terms proportional to ^ are numerically much smaller than those 
proportional to 4. This is not surprising, since the e.m. coupling of the photon to the 
quarks corresponds to a magnetic transition, and therefore it is inversely proportional to 
the quark mass. 
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A final remark concerns the photon energy. In eq.(|H|) we have allowed k° to vary 
in the range [0, Mb/2]; however the integral diverges at the lower limit k° = 0. This 
result is unphysical since it would correspond to a zero energy photon in the final state. 
Formally, this divergence would be canceled by radiative corrections to the formulae fll5D 
and (0). On the other hand, one should take into account that at future experiments, 
e.g. at the SLAC 5-factory, the smallest measurable photon energy is of the order of 
50 MeV; therefore, it is a reasonable assumption to cut off from the integral the small 
photon energies, and, at the same time, to neglect radiative corrections. 

In our calculation, the effect of the unphysical divergence begins around a photon 
energy k° ~ 350 MeV, and therefore we use this value as a lower bound for the photon 
energy. 

In fig. 3 we plot the photon spectrum for the decay (H); the differential distribution has 
a peak around 1.5 GeV, which should render it quite accessible to experimental analyses. 

For the decay width we obtain the result T(B~ — > ^"iyy) = 3.7 10~ 19 ( 3x v |" ''_3 ) 2 GeV, 
which corresponds to 

B(B~ -> n~v^) = 0.9 10~ 6 . (23) 

This result is obtained with a cut-off A = 350 MeV in the photon energy. As it can 
be seen from fig. 3, the uncertainty related to this choice should not be significant. For 
example, putting A = 100 MeV the result for the branching ratio increases by less than 
10%. 

The conclusion we can draw from this result is that the relativistic quark model gives 
predictions in agreement with the expectations discussed in ref. ||]; the weak radiative 
decay B~ — > pTv^ has an appreciable rate and might be observed in the near future. 

As for on the theoretical uncertainties of the result (0), a part from the energy cut-off 
they mainly come from the choice of the B meson wave function. A possible estimate of 
this theoretical error consists in assuming a different wave function. A commonly used 
quark momentum distribution inside the B meson is given by the ACCMM model ||12|| ; 
in our notations it corresponds to the gaussian wave function: 

u B (k l )=2fk l \(^y'\ Xp (-§,). (24) 

The parameter Pp is related to the heavy quark average square momentum: < p 2 >= f-Pf . 
In a recent analysis [H| this parameter has been fitted using experimental data for the 
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inclusive decay B — > X c £u£ fTJj] with the result: Pp = 0.51 GeV^. Using (B3) with such a 



value of Pp in the previous formulas, one would obtain B(B — ► // 1^7) — 0.8 ■ 10 . This 



result suggests that the estimate given in eq. (p3|) is rather accurate; the rate for the weak 
radiative decay E>~ — > is large enough for a measurement at future accelerators. 
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Figure Captions 



Figure 1 

Diagrams describing the decay B~ — > fj,~Vf/y. 
Figure 2 

The wave function us(ki) as obtained by the QCD relativistic quark model. 
Figure 3 

Predicted photon energy spectrum. 
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